Phase locking dynamics of dipolarly coupled vortex-based spin transfer oscillators 
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Phase locking dynamics of dipolarly coupled vortices excited by spin-polarized current in two 
identical nanopillars is studied as a function of the interpillar distance L. Numerical study and 
analytical model have proved the remarkable efficiency of magneto-static interaction to achieve 
phase locking. Investigating the dynamics in the transient regime towards phase locking, we extract 
the evolution of the locking time r, the coupling strength /x and the interaction energy W. Finally, 
we compare this coupling energy with the one obtained by simple model. 



Injecting a spin-polarized current through magnetic 
multilayers leads to new interesting physical phenomena 
named spin-transfer effect. These interactions between 
the spins of charge carriers and local magnetic moments 
create an additional torque exerted on the magnetization 
As a result, a complex spin-transfer-driven magnetic 
dynamics comes out with characteristic bifurcations of 
the Poincare-Andronov-Hopf type and limit cycles arises 
in this highly non-equilibrium medium. The diversity of 
these new effects is especially true for systems of interact- 
ing nanomagnets, penetrated by spin-polarized current. 
One of the novel effect is the current-driven magnetiza- 
tion oscillations that might lead to tantalizing pos- 
sibilities for new nanoscale microwave devices with the 
frequencies tunable over a wide range using applied cur- 
rents and fields. While many crucial advances have been 
made in the fabrication and understanding of such Spin 
Transfer Nano- Oscillators (STNO), there remain several 
critical problems yet to be resolved, in particular, the low 
microwave power and quality factor of a single STNO. 

To tackle these issues, particular attention has been re- 
cently focused on vortex STNOs that could present a sig- 
nificant output power [3] , a very small spectral linewidth 
[1] and/or large frequency agilities at zero field [5]. More- 
over, several encouraging experiments have been reported 
on the vortices phase-locking through exchange interac- 
tion \m and synchronization to external microwave cur- 
rent [7|. Beyond these practical interests, a magnetic vor- 
tex and its dynamical modes [8], notably the gyrotropic 
motion of the vortex core, is a model system to investi- 
gate deeply the physics of the spin transfer torque acting 
on a highly non uniform magnetic configuration [l^ . 
Collective gyrotropic modes is a mean to improve drasti- 
cally the spectral coherence of any oscillator system . 
Similarly, vortex based systems can be chosen as a new 
playground to investigate the influence of the magneto- 
static interactions on vortices collective behaviour. 

Magnetostatically coupled vortices collective dynam- 
ics have been studied both experimentally and theoret- 
ically for the case of low amplitude oscillations excited 
by means of external RF magnetic field [12-16] and spin- 



amplitude steady oscillations. The fundamental reason 
for that is the hypothesis of low amplitude oscillations 
near the centers of nanodots used by these models. A 
strong consequence of this approximation is that mathe- 
matically it allows to use of the ratio of the vortex orbit to 
the disk radii as a small parameter. However in the case 
of the large-amplitute oscillations phase-locking, such lin- 
earization is due neither for the of vortex STNOs, nor for 
the uniform ones [I^] • In the letter we propose an original 
model for the coupled vortices dynamics without using 
this assumption. This model provides an expression for 
the coupling energy with the parameters of the transient 
process, which can be directly determined either through 
micromagnetic simulations or by experiment. 

The studied system is made of two identical nanopillars 
with diameters 200 nm, each of them being composed by 
a free magnetic layer, a non magnetic spacer, and a since 
synthetic antiferromagnet (SAF) polarizer which gener- 
ates a perpendicular spin polarization pz. In our simula- 
tions we consider these layers only by the value of spin- 
polarization like in [17], since SAF polarizers that are 




polarized current [17|, ll8| . However all these models are 
not applicable to the case of interest, i.e. for the large 



FIG. 1: (Color online) Schematic representation of two inter- 
acting spin-transfer oscillators. Each pillar is composed by a 
free magnetic layers with vortex, a non magnetic spacer, and 
a SAF polarizer. Red arrows indicates the direction of spin 
polarization created by the polarizer. The nanopillars have 
a diameter D = 200 nm and are separated by a distance L. 
The parameters Xi and X2 define the cores positions. 



widely used in vortex STNO experiments, have a neg- 
ligible magnetostatic field, and thus practically do not 
influence on the vortices dynamics. 

A free layer is = 10 nm thick NigiFeig and has 
a magnetic vortex as a ground state. The magnetic 
parameters of the free layer are: the magnetization 
Ms = 800 emu/cm^, the exchange energy A = 1.3 x 
10~^ erg/cm and the damping parameter a = 0.01. In 
order to be above the critical current, a spin polarization 
P of 0.2 and a current density J of 7 x 10^ A/cm^ have 
been chosen. The initial magnetic configuration is two 
centred vortices with same core polarities and same chi- 
ralities. The micromagnetic simulations are performed 
by numerical integration of the LLG equation using our 
micromagnetic code SpinPM based on the forth order 
Runge-Kutta method with an adaptive time-step control 
for the time integration and a mesh size 2.5 x 2.5 nm^. 

In this work the evolution of the phase-locking dynam- 
ics as a function of the interpillar distance has been stud- 
ied. Therefore a series of micromagnetic simulations with 
different distances L (50, 100, 200 and 500 nm) has been 
performed. The results of the simulations are then ana- 
lyzed to extract the radius of the vortex core trajectory in 
each free layer as well as the phase difference tjj between 
core radius- vectors as a function of time. 

In Fig. [21 the simulations results for L = 50 nm are 
presented. The vortices transient dynamics can be di- 
vided in two regimes. At t = the spin torque is switched 
on and thus the radii of both cores trajectories increase 
towards their equilibrium orbits for about 300 ns (see 
Fig. [2^). The phase difference between the two radius 
vectors shown in Fig. [JJd remains constant and equal to 
— TT because of repulsive core-core interaction. The sec- 
ond regime begins when the two cores have reached orbits 
close to their steady ones. From this stage, both the in- 
tercore distance and the phase difference (see Fig. [2^ and 
b) show large oscillations indicating the beginning of the 
phase locking. 

The second regime is one of main interest for this work 
since the couphng energy can be extracted using the anal- 
ysis of the cores motion in this transient regime (indi- 
cated by the square in Fig. [JJd). During this time range, 
the phase difference ip can be identified as being a low 
frequency damped oscillation described by the following 
expression: 



sin(m + C2). 



(1) 



As shown in Fig. [2t, the fitting is done for the time 
window between 500 and 800 ns in which the mean or- 
bit radii have reached the common equilibrium value Xq . 
From the fitting, one can extract for L = 50 nm a fre- 
quency Q equals to 40.134 MHz and a phase locking time 
of 82.59 ns. 

The parameters extracted from the fitting procedure 
for all the interdot distances are summarized in table HI 
One needs to note that the phase locked equilibrium orbit 
radius Xq does not vary much with L. 

To derive the coupling energy between the oscillators 
from the simulations, we have developed an analytical 
model based on Thiele equations [20| coupled through the 
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FIG. 2: (Color online) Micromagnetic simulations for L = 50 
nm of the phase locking dynamics. Evolution as a function 
of time t of the vortex core orbital positions Xi and X2 (a) 
and the phase difference ip (b). In (c), a zoom of the phase 
difference is presented for the time window in which the 
fitting with Eq. ([T]) has been done. 
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FIG. 3: (Color online) Phase difference '0 as a function of 
time t for different interpillar distances L = 100 (a), 200 (b) 
and 500 nm (c). 



L (nm) 


n (MHz) 


T (ns) 


^0 (nm) 


50 


40.134 


82.59 


63.59 


100 


28.305 


85.28 


62.46 


200 


17.183 


89.57 


61.82 


500 


7.018 


90.13 


61.53 



TABLE I: Values of the phase difference frequency Q, the 
phase locking time r, and the steady- state radius of the core 
motion Xq at different interpillar distances L using the ex- 
pression of Eq. ([Tj) 



dipolar interaction energy Wint 2ll-[23|. Due to the sys- 
tem symmetry, the interaction energy can be expressed 
as Wint = ciiXiX2-\-biyiy2 that can be reformulated using 
the cores positions Xi and X2 as (2^ : 

Wint(Xi,X2) = /iiXi • X2 + /i2(^i3:2 -ym)- (2) 

where /ii^2 are the interaction parameters and xi^2, yi,2 
are core coordinates. The second term of Eq.(j2]) is ne- 
glected in our study, since it corresponds only to fast 
oscillations at double frequency of the gyrotropic modes 
and thus is averaged over the low frequency dynamics 



which is responsible for the phase locking. As a conse- 
quence, the expression for the interaction energy can be 
written as Wint 0^1,^2) = • X2. 

The two Thiele equations of the core dynamics consid- 
ering both the spin transfer torque and the interaction 
between the two oscillators are: 



G 



(e, X Xi,2) - fei,2(X)X - Z)Xi,2 - FsTT 



■ int 



= 

(3) 

The first three terms are the conventional forces: the 
gyrotropic force with Ge^ = —27Tp^^ez^ the confining 

force with k{X.) = uoG (^1 + where the gy- 

rotropic frequency is ujq = ^jMsh/R and the dyadic 
damping D = arjG^ 7^ = ^ In (^^^ + |, where Iq = 

2^^2 ' The fourth term Fstt is the spin transfer force. 

For the case of uniform perpendicularly magnetized p< 
larizer Fstt = 7rjajMsh{ez x X) = ^(e^ x X) 
where the spin torque amplitude is aj = PJ/Mgh with P 
the spin polarization and J the current density. The last 
term describes the interaction forces and is expressed by 

F.nt(Xi,2) = -SI^,nt(Xi,X2)/SXi,2 = -/iX2,l. ThcSC 

Thiele equations can be reformulated in polar coordinates 
as (using ar] <C 1): 



Xi 

X2 
X2 



arikiiXi) -\- a X2 , . , 

) +g^(sinV^-«rycosV^) (4) 

G / - ^^(smV^ + arycosV^) (5) 



f (l/{rar]) - ^/l/{Tar])^ - n{L)^^ Xl In Fig. H the 
evolution of this interaction energy Wint with the inter- 
pillar separation distance D12 derived from micromag- 
netic modelling is displayed by blue square dots. The 
best fit we obtain if for an energy decay law to be 
Di2 '^. In comparison, in the case of small core ampli- 
tudes m, [l3| this decay law has been found as D^2^ 
however without being confirmed experimentally ^7|] . 
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FIG. 4: (Color online) Absolute values of interaction en- 
ergy Wint as a function of interpillar separation distance 
D12 — 2R + L obtained from the micromagnetic simulations 
(blue square dots) and from the macro-dipole model (red line). 
Inset: Evolution of the phase locking rate l/r versus D12 
(purple line). 



ip — auo 



xl 



Xl X2\ 1^ , fX2 Xl 

(6) 

These two equations of core motion and the equation of 
phase difference provide a complete dynamical descrip- 
tion of the phase locking. By linearizing Eq. ()4|5|6p with 



the assumptions that <C 1 and e 
obtain: 



< 1, we 



X1+X2 



e = 2ar]{il — (jjQar1)e ^ jlil) 
= 4:{fi — ujoarl)e — 2ar]fitlj 



(7) 
(8) 



where we used ji = (i/G and ro = Xq/R. The two equa- 
tions ([7j) and dHI) are linear and their eigenvalues are 

Ai,2 = ariujQarl±^ a^rfujl^arlY ^ - Afiuj^arl. (9) 

consider the case of periodic solu- 

< jl < 



First we 
tions when 



I (^uJoarQ -\- uJoarQ^/T^^o?rf^ . The eigenvalues can be 



thus written as: 



l/r = ar]ujQarQ 

o2 2 2 2/ 2\2 



4/i^ + AfiuJoarQ 



(10) 

(11) 



The important result of this study is that Eq. ([TC 
and (pT]) allows to connect the coupling parameter /i 
with the phase locking parameters, i.e. Q and r, 
obtained through micromagnetic simulations. Conse- 
quently, an expression of the time-averaged interac- 
tion energy Wint takes the form: Wint{L) = /iXg = 



To get more insights about the origin of this large 
coupling interaction, the values of the interaction energy 
Wint obtained by simulations are compared with the ones 
derived from a simple model of two interacting macro 
dipoles, concentrated at the dot centres and rotating at 
a frequency CjOq. In such a case, the magnetic dipole in- 
teraction energy W^^^ is defined as 



(Ml •M2) 



^12 



3(Mi-Di2)(M2-Di2) 



^12 



(12) 

where D12 = 2R-\-D. The in-plane magnetization Mi^2 is 
perpendicular to the radius-vector of core position, thus 
one can write Mi^2 = C(^i,2 x e^), where ^ is a constant, 
that has been numerically calculated: ^ ^ 5.6 G/nm. 
The interaction energy can be rewritten in the following 
form: 



with A 



AX1X2 

low frequency oscillations 



2I)?o 



B 



2 I^?o 



- 5X1X2 cos((/?i + (P2) 

^ V ' 

high frequency dynamics 

(13) 

As far as the phase 



locking dynamics is concerned, the second term in (p!3|) 
corresponding to high frequency oscillations is averaged 
to zero and thus one can express the mean interaction 



energy W^^^ ^ in the macro-dipole approximation: 



int 



X1X2 



m — d 



X1X2 



(14) 



In Fig. |4]we observe that for small interpillar distances 
Wint differs significantly from the macro dipole energy 



W^^^~ , this difference demonstrates the importance of 
the magnetic quadrupole and higher multipoles for the 
phase- locking dynamics. 

Coming back on the Eq. ^ and dU, 

a second regime has to be considered when 
jl < \ {oJoarl - ujQarl^l - a^rf^ or ft > 

\ [ujQarQ + ujQavQ \/l — a^r]'^^ . In this case, the so- 
lutions are aperiodic oscillations and it strongly impact 
the main features of the phase locking, notably the 
phase locking rate 1/r. Indeed, in the regime of periodic 
oscillations, this phase locking parameter is almost 
independent on the coupling strength /i for interpillar 
separation distance D12 values as large as 1600 nm (see 
inset of Fig. |4]). We emphasize that the weak variation 
of the phase locking rate obtained in the micromagnetic 
simulations (see values in table H]) is in fact solely due 
to the small variations of the steady orbit radii Xq with 
the interpillar distance L as expected from (p!Q|) . On the 
contrary, in the aperiodic regime, the phase locking rate 
l/r depends strongly on the coupling strength /i with 
the following expression: 



l/r = rjujoarl - ^J r]'^uj'^{ar'^y + 4/i2 - AfiuJoar'^ (15) 

Using the value of the coupling parameter jn that can 
be extracted for very large interpillar distance L through 
the macro-dipole approximation, we obtain a very rapid 
decrease of l/r with interpillar distance L and eventu- 
ally a phase locking time that tends to r — > 00 for large 
distances. It is important to note that interaction en- 
ergy Wint becomes of the same order of magnitude as 
the room temperature thermal energy kT at the inter- 
pillar distances L of about single STNO diameter, thus 
the role of thermal effects in the phase-locking of vortex 
STNOs has to be properly investigated. 

In conclusion, we have demonstraded an efficient 
phase-locking between two STNOs through dipolar 
mechanism. We have succeeded to provide an accurate 
expression of the interaction energy between two vortices 
based STNOs by comparing micromagnetic simulations 
to predictions of an analytical model based on coupled 
Thiele equations with dipole-dipole interacting forces. A 
major result is that the phase locking time r is almost 
independent on the separation distances for values up to 
1.6 /im before it increases very rapidly at larger distances. 
We emphasize also the critical importance of higher order 
multipole terms for a correct description of the interac- 
tion energy, especially at shorter separation distances. 
Finally, our investigation opens ways to design some op- 
timized STNO ensembles for synchronization which is a 
crucial step toward the development of a new generation 
of RF devices for telecommunication applications. 
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